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Abstract 



Relations between Wigner functions and the source functions used in 
t models of Bose-Einstein correlations in multiple particle production are 

derived and discussed. These relations are model dependent. In particular 
it is important whether the particles are emitted simultaneously and if not, 
whether the production amplitudes corresponding to different moments of 
1 time can interfere with each other. 

a: 

^ 1 Introduction 

jz - 

When discussing Bose-Einstein correlations among identical pions 1 produced 
in multiple particle production processes, one often uses the source function 
. S(X, K) (cf. e.g. PP and references quoted there) related to the single particle 

density matrix in the momentum representation by the formula 2 

p(p,p') = J d*Xe^ x S(X,K), (1) 

where q = p — p' and K = (p + p')/2 are four- vectors constructed from the 
on-shell single particle four momenta p and p . The four components of X are 
integration variables and, therefore, there is much freedom in their interpreta- 
tion. It is usual, however, to identify S(X, K) with the space-time distribution 

'Partially supported by the KBN grant 2P03B09322 

1 Here and in the following we call for definiteness the identical bosons pions the results, 
however, are valid for any identical bosons 

2 Often an equivalent formula with the density matrix replaced by the average 
< a^iipi >= p(p', p) is used. Also some authors include on the left hand side a factor VEE' . 
Since this factor does not affect our argument and is easy to introduce at any stage, we skip 
it. 
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of the sources producing the pions with momentum K and then X is interpreted 
as the space-time position of the source. A source functions combines conve- 
niently the information and/or prejudice about the space-time (X) distribution 
of the sources and about the momentum distribution of the produced particles. 
It yields the density matrix, which can be used to get the (two- as well as more 
particle) correlation function. Comparing the predicted correlation functions 
with experiment one can improve the model used to find S(X,K), fix its pa- 
rameters etc. Note that the density matrix on the left hand side of formula (JTJ 
does not depend on time. Since after freezeout the particles propagate freely (we 
do not discuss the final state interactions here), this means that for times after 
the freeze out period the elements of the density matrix should be interpreted as 
the matrix elements of the time dependent (Schrodinger picture) density opera- 
tor between the time dependent states |p,i) = e <B ^*|p), or equivalently as the 
matrix clement of the time independent (interaction picture) density operator 
between the time independent states |p). 

Formula JJJ looks similar to the formula relating the Wigner function W(X, K) 
(cf. e.g. [2] and references contained there) to the density matrix: 

P(P, p', t) = J d 3 Xe-^ X W(X, K, t). (2) 

In fact the building block of the source function was originally a Wigner function 
(cf. formula (7)). Because of the similarity between formulae and J5J, 
the source function is often referred to as a Wigner function, a kind of Wigner 
function, a pseudo- Wigner function etc. (cf. e.g. ref. 4 , from which all 
these names have been taken). A discussion of the actual relation between the 
source function and the Wigner functions, however, seems to be missing in the 
literature. Let us begin with some general remarks. 

Formula J5J) contains the time dependent density matrix, i.e. the time de- 
pendent density operator in the representation of the time independent states 
|p) . Thus 

p(p,p',i)=e-^V(P,P')- (3) 
Relation follows from (J2J), if we put 

S(X,K) = W{X,K,t )6(X -t ). (4) 

This choice can be always made and has a simple physical interpretation. Choose 
a time to after freezeout, when all the pions are already present, but none have 
been registered yet. The state at to yields the initial condition from which any 
distribution at registration can be calculated. Assuming that all the pions ap- 
peared simultaneously at to may be poor physics, but this does not change the 
fact that the predictions at registration time are correct. There is an infinite 
choice of other source functions, which also satisfy relation QJ. Since this rela- 
tion is the only link between a source function and physical reality, there can 
be no objection of principle against the choice Q. From the practical point 
of view, however, another source function, which can be calculated from some 
model, may be more convenient and, of course, it is just as good, if it reproduces 
the density matrix equally well. Note that the identity of functions in J3J does 
not mean that their argument have the same physical interpretation. In fact, 
on the left hand side X and K denote the position of the source in space-time 
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and the four-momentum of the pion, while on the right hand side X = |(x+x') 
and K = |(p + p') are defined by reference to the arguments of the density 
matrices of the pion in the coordinate representation and in the momentum rep- 
resentation respectively, while to is the time. It is often reasonable to assume 
that X and K on both sides are good approximations to the position at t = to 
and momentum for t > to of the pion, but this is certainly not always the case. 

Let us conclude this section with some more remarks. If one assumes that all 
the particles were created simultaneously at some common freezeout time t = to, 
then a "realistic" source function should be proportional to S(Xq — to) and the 
proportionality coefficient is the Wigner function as in |J2J| . The source function 
cannot be equal to a Wigner function, because the dimensions are different. 
The delta function in formula brings into the dimension the necessary factor 
1/time. What could be equal to the source function is what we will refer to as 
the differential Wigner function defined by the relation 

W(X,K,t ) = W dto pt,K,t )/dt . (5) 

Here Wdt is the Wigner function at time t > to for the particles from sources 
created in the time interval dto around time to- Note that we are using the 
Wigner functions normalized to the numbers of particles and not to one. When 
the amplitudes of the pions produced at different times add incoherently, formula 
iPJ can be rewritten in terms of the differential Wigner function W as: 

p(p,p') = J d 4 Xe^ x W(X,K), (6) 

where the fourth component of the vector X = (X, Ao) is to- Comparing this 
formula with formula we find a particular solution for the source function 

S(X,K) = W(X,K). (7) 

We stress that this is only one out of the infinity of source functions, which 
when substituted into formula Q yield the correct density matrix and that for 
the sources incoherence in the creation time has been assumed. 



2 GGLP models 

In the seminal paper of G. Goldhaber, S. Goldhaber, W. Lee and A. Pais 5 the 
assumption was made that pion production amplitudes at different space points 
are incoherent. In this simple case there is no need to distinguish between the 
sources and the pions at t = to- A natural generalization (cf. e.g. [S]) was to 
assume more generally that pions production amplitudes at different space-time 
points are incoherent. Further we refer to such models as GGLP models. 

Let us begin with the simpler case, when all the pions appear simultaneously 
at t = to- For t > to they are assumed to propagate freely. The time independent 
density matrix for t > to is: 

P(P,P')= / d 3 x{p,t\x,t;t )p{x,t ){x,t;to\p't). (8) 

Here |x, t;io) = exp[— i(t — to)Ho\\x) is the time dependent state, evolving ac- 
cording to the free particle Hamiltonian Ho from the initial state |x), which 
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corresponds to a pion localized at x at time to- One could also say that in 
space-time the position of the source is (x, to)- Let us note for further use the 
formula 



(p,t\x,t;t ) = e iE W t °{p\x), (9) 

where E(p) = \/p 2 + rn%. Formula (JSJ) is perfectly acceptable from the point 
of view of quantum mechanics, in spite of the fact that according to the picture 
described in ref. [5] a pion with momentum p is produced at point x, which 
would contradict the uncertainty principle. Since the density matrix is diagonal 
in x (in the x, which labels the incoherent states, not in position at time t) and 
t = t' = to , we can identify the components of the vector x with Xi , X2, X3 
and the parameter to with Xq. The matrix element can be rewritten as 

f d 3 X 

P(p,p') = J j^exp(iqX)p(X). (10) 
Using © and comparing with this yields 



W(X,K)=>£-L, (11) 



P(X) 
(2tt) s 

The source function is given by formula Q, because by assumption all the 
particle production happens at t — to- The components of X are the three 
components of the position vector of the point where the particle was produced 
and its time of production (freezeout time) . Since after freezeout the momentum 
vector of each particle is conserved separately, the momenta entering K can be 
taken at any time not earlier than to- The actual time of observation does not 
occur in the formula at all. 

Note that for this class of models the Wigner function does not depend on 
momenta. With the present interpretation this would lead to disaster when 
comparing with experiment. The Goldhabers Lee and Pais |S1 proceeded dif- 
ferently. They were interested in n-pion exclusive states for n = 4, 5, 6 with 
no more than two pions of a given charge. Using matrix (jSJ) they constructed 
the probability distributions P(j>i, ■ ■ ■ ,p n )- These were interpreted as momen- 
tum distributions only after the sets of momenta contradicting energy and/or 
momentum conservation were projected out. In practice they interpreted the 
P-s as integrands of the corresponding phase space integrals. This could not 
be generalized to inclusive processes, led to rather complicated integrations and 
gave results, which at that time were considered acceptable, but which soon 
were shown to be in violent disagreement with experiment |7j. Most of these 
difficulties were overcome by Kopylov, Podgoretsky and coworkers as described 
in the next section. 

Let us consider now the more general case, when the pions are produced at 
different moments of time, but there is no interference between the amplitudes 
corresponding to production at different moments. Then at times larger than 
the largest freezeout time the density matrix is given by an integral over the 
freezeout time to of the (differential version of the) density matrix I] 10(1 corre- 
sponding to particles produced at given time to'. 

P(P,P') = J '0-exp(iqX)p(X). (12) 
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Note that in order to keep the dimensions right p(X), which is a density in 
space, has been replaced by p(X) which is a density in space time. In fact, 
p(X) are diagonal elements of the time independent single pion density matrix 
in the |x, t;to) representation, while p(X) are the diagonal elements of the cor- 
responding differential density matrix related to p(X) by a formula analogous 
to JSJ. Comparing with (JTJ it is seen that one solution for S(X,K) is 



Since there is incoherence in time, the relation between the source function and 
Wigner's function is S(X,K) = W(X,K) 0. 

This seems to be the perfect situation: the source function reproduces ex- 
actly the distribution of (incoherent) sources in space-time and coincides with 
the differential Wigner function. However, l|13(l is only one of the infinity of 
source functions, which yield the same density matrix. A quite different one 
can be obtained e.g. by calculating the Wigner function corresponding to the 
density matrix p(p 1; p 2 , t) and substituting it into formula (@J. We can conclude 
only that for this class of (unrealistic) models there is an acceptable source func- 
tion, which can be interpreted as a differential Wigner function. Therefore, it 
can be argued that this particular source function is easier to find than others. 

3 Including coherence in space 

An important extension of the GGLP approach is to assume that at any given 
creation time there may be coherence between the amplitudes of the pions pro- 
duced at different points in space, though pions produced at different times to do 
not interfere. Models of this type, also in the more general version with coher- 
ence in time, have been introduced by Kopylov and Podgoretsky 8 a . In terms of 
sources this may be quite complicated. For instance, in the model described by 
Kopylov and Podgoretsky |Hj all the sources are produced simultaneously and 
live for some time. The averaging over the energy of each source (not of the 
pion!), however, kills the interference between the amplitudes for the production 
of the pion by the same source at different moments of time. Therefore one can 
replace each original source by a set of sources distributed in time and produc- 
ing the pions instantly. These are not quite satisfactory models, because it is 
difficult to explain why the creation amplitudes from neighboring time moments 
should not interfere, while those from neighboring space points do. Moreover, 
relativistically the creation processes which are at different times in one Lorentz 
frame may be simultaneous in another Lorentz frame. As seen from formula 
fllfl. however, even in the simplest case of simultaneous production, such models 
can reproduce correctly any single particle density matrix and consequently any 
single particle momentum distribution without doing phase space integrations. 
Thus they can be applied to describe inclusive processes. 

Including interference in space and summing without interference over time 
corresponds to the following generalization of formula JHJ) 



p(p,p') = dt d 3 x / dV(p,i|x,i;i }p(x,x',i )(x',i;to|p',i}. (14) 



8{X,K) 



~P{X) 



(13) 
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There are many equivalent ways of rewriting this relation. For instance, one 
may use a representation \a), where the (differential) density matrix p(x, x',to) 
is diagonal. Then l|l 4|) gets replaced by 

p(p,p') = J dt J dae I9ot0 (p|a}p(a;t )(a|p'}. (15) 

It is instructive to compare this model, to the GGLP models. At to the localized 
states |x) = \x,t ;t ) have been replaced by the states |a) = \a,t ;t ). The 
plane waves (p|x) have been replaced by the functions (p|a), which are known 
under a variety of names: as waved packets, as sources, or as currents. The 
source function for such models was introduced by Pratt 3 . His statement that 3 
" S(X, K) can be identified as the probability of emitting a pion of momentum 
K from space-time point X" is, however, only an approximation. 

Kopylov and Podgoretsky 8 assumed that (a, to) is a point in space-time, 
so that formula (|15fl can be rewritten as 

p(p,p') = J d 4 x (p^ xo )p(xo)^ xo \p')e^ to . (16) 

This approach yields an alternative method of describing the distribution of 
pions in phase space. In the integrand (p\ip Xo ) is the probability amplitude 
for finding a pion with momentum p produced by a source labelled xo- When 
the states \ip Xa ) correspond to particle well localized in space, p and Xo give a 
reasonably good description of the position of the pion in phase space. One could 
object that this is only a rough description, but the same is true for the Wigner 
function: K and X give only approximately the momentum and position of the 
pion. An exact determination of the pion position in phase space is possible 
only in classical physics. 

In ref. [S] the states \a) = \tp xo ) were supposed to be related by space time 
translations so that 



ip Xo (x) = (x\tp Xo ) =ip(x-x ), (17) 
(p|^ ) = e**V(p), <£(P) = (Pl^o). (18) 

The second formula is, of course, equivalent to the first. Another choice has been 
made in the "covariant current ensemble formalism" ^U], [TT], . There 

each source is labelled by a position in space-time Xq and a four- momentum po- 
Usually xo and p denote the centers of the wave packet in ordinary space and 
in momentum space respectively. The time component of xq is to - the freeze 
out time of the wave packet. The time component of po is calculated from the 
condition p 2 , = m 2 , where m can but does not have to 0] be the pion 
mass. The nice feature of this parametrization of the wave packets is that one 
can substitute a classical trajectory Xo(i),p (i) for the source and remain in 
agreement with the Heisenberg uncertainty principle for the pions. The choice 
for the scalar products is 

<PlVWo>=e^°j(— )• (19) 

TO 

3 This has been changed to our notation. Pratt has written g,p,x where we have written 
S,K,X 
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Thus, the sources are related by space-time translations and when each of the 
currents j is considered in its rest frame where po — (m,0), they are identi- 
cal. Formulae (|17() and l|19l) have been applied also to relativistically covariant 
models cf. e.g. Pj. 

Another way of rewriting relation (|14|) is 

P(P,P') = J d 4 Xexp(iqX) J ^-^3 exp(-iK-y)p(x,x / ,to), (20) 

where X = (x + x')/2, y = x — x' and Xq — to- According to formula a 
possible choice of the source function is 

S(X, K) = J -j^L expHK • y)p(x, x', to). (21) 

The differential Wigner function is the Wigner transform of the differential 
density matrix as it should. The relation of this particular source function to 
the Wigner function is again given by formula JJJ). We conclude that among 
the infinitely many source functions which give the same density matrix in the 
momentum representation there is one, which can be related to the Wigner 
function as described. Note that this source function depends on K, but does 
not depend on Kq. 



4 Including coherence in space and time 

Finally let us consider the case, when neither coherence in space, nor coherence 
in time is assumed. Then the density matrix is 

p(P,p') = / d*x J d i x'(p,t\x,t;t )p{x,x')(x.',t;t' \p , ,t), (22) 

where x = (x, to) and x' = (x', t' ). For an application of a formula of this type 
see e.g. ref. [0], where ^2 F Mf(x)Mjp(x') stand for our p(x,x'). 
Formula (JUJ) can be rewritten as 

p(p,p') - Jd 4 X J d 4 yexp^X + iKy]p(x,x') /(2ir)\ (23) 

where X = ~(a; + x') and y = x — x'. Comparison with formula gives as one 
of the solutions for the source function 

S{X,K) = J Ay^ v [iKy]p{x,x'). (24) 

The right hand side is rather remote from what one usually calls a Wigner 
function. The differential Wigner function was at least proportional to a Wigner 
function, though only for the particles from sources produced in the time interval 
dto around the time t — tg. Function (|24|l is proportional to the contribution 
to the Wigner function from the interference of the production amplitude in 
the time interval dto around t = to with the production amplitude in the time 
interval dt' around t = t' , integrated over t' — to at fixed t' + to- 
Formula l|22|) rewritten in terms of wave packets reads 
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p(p,p') = J dto I dt' I da{p,t\a,t;t )p(a,to,t' )(a,t;t' \p',t). (25) 

Note that the state \a,t;to) may, but does not have to, be connected to the 
state \a,t;t' ) by a smooth hamiltonian evolution. The source at time to can 
be, as well, something quite different from the source at time t' . On the other 
hand, the evolution of \a, t; to) in the time t for t later than the freezeout of the 
particle is the ordinary free particle evolution. The various models are defined 
by the choice of a, of the states \a, t; to) and of the weight function p(a, to, t' ). 
Examples can be found in refs. [S], 0]. 



5 Conclusions 

A source function cannot be equal to a Wigner function, because they have dif- 
ferent dimensions. Moreover, for a given state of the system its Wigner function 
is well defined, while its source function is not. The problem is, therefore, to 
chose some special source function and try to relate it somehow to a Wigner 
function. 

When all the particles are produced simultaneously and when this is assumed 
to mean that the source function is proportional to a delta function in time, the 
proportionality coefficient is unambiguously defined as the Wigner function of 
the pions at the time of freeze out (J2J . When it is assumed that the production 
amplitudes at different times add incoherently, one can use the source function 
proportional to the differential Wigner function as given in (|13fl . When the 
production amplitudes at different times interfere, a source function can be 
related to a piece of the Wigner function as given by formula l|24|l and explained 
below this formula. In this case, however, the use of wave packets (or sources, 
or currents) may be a more natural way to analyze the phase space distribution 
of pions. 

If, as is often done, a model is defined by postulating the source function 
(cf. for examples), the question about the relation of this source function to 
a Wigner function cannot be answered without making an assumption about 
the coherence or incoherence of the amplitudes for pion production at different 
times. When the source function is proportional to a delta function in time, one 
can relate it to a Wigner function by relation When it does not depend on 
Ko, one can use formula Q. In the later case there is no guarantee that the 
production process was such as the differential Wigner functions suggest. 
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